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Abstract. Under the assumption that jets in active galactic nuclei are powered by accretion and the spin of the 
central supermassive black hole, we are able to reproduce the radio luminosity functions of high- and low-excitation 
galaxies. High-excitation galaxies are explained as high-accretion rate but very low spin objects, while low-excitation 
galaxies have low accretion rates and bimodal spin distributions, with approximately half of the population having 
maximal spins. At higher redshifts (z ~1), the prevalence of high accretion rate objects means the typical spin was 
lower, while in the present day Universe is dominated by low accretion rate objects, with bimodal spin distributions. 



1. Introduction 

Active galactic nuclei (AGN) are galaxies with signs of 
non-stellar activity in their centres, believed to be pow- 
ered by supermassive black holes (SMBHs). AGN produce 
jets which are observable at radio frequencies. For a given 
accretion rate, AGN can produce jets which vary in radio 
luminosity by several orders of magnitude. Conversely, for 
a given radio luminosity, some AGN are found to have 
very high accretion rate, while others are found to have 
very low ones. 

The hidden variable behind this 'radio loudness' of 
AGN has often been connected with the spin of the SMBH, 
a, and this is known as the 'spin paradigm' (e.g. Wilson & 
Colbert 1995; Sikora et al. 2007). Indeed theory and sim- 
ulations suggest that spinning SMBHs should be able to 
power considerable jets (e.g. Blandford & Znajek, 1977, 
Hawley & Krolik 2006, Tchekhovskoy et al. 2010), and 
indeed observations of the most powerful radio galaxies 
require an extremely efficient mechanism for the produc- 
tion of jets (e.g. Punsly 2007, Fernandes et al. 2011). 

Recently, Fender et al. (2010) have noted a lack of cor- 
relation between the reported measurements of jet pow- 
ers and black hole spin for galactic black holes in X- 
ray binaries. As we discuss in detail in Martmez-Sansigre 
& Rawlings (2011, hereafter MSR11), the uncertainties 
in both measurements are very large, so that the lack 
of observed correlation does not provide robust evidence 
against the spin paradigm (see Section 7.2 of MSR11 for 
a detailed discussion). 

In this work we infer the cosmic distribution of SMBHs 
spins under the assumption that spin is the major variable 
explaining the huge variations in jet power. Indeed mod- 
els and simulations of spinning black holes can provide 



enough variation in jet power to explain most of the vari- 
ation in radio loudness of AGN (see Section 3 of MSR11). 

Under this assumption, we derive the spin distribu- 
tions for SMBHs with high and low accretion rates, and 
infer the cosmic spin history. Throughout our work we 
concentrate on the most massive black holes, with m»> 
10 s M Q . This is achieved by working only with radio lumi- 
nosity densities > 10 23 W Hz -1 sr _1 , where all the AGN 
black hole masses are known to be > 10 8 M Q (McLure et 
al. 2004, Smolcic et al. 2009). 

1.1. Assumptions 

The bolometric power available for radiation is given by: 
iboi = e(a)m.c 2 (1) 

where £t>oi is the bolometric luminosity due to radiation, 
m. is the rate of accretion of mass onto the SMBH, and 
the term e(a) is the radiative efficiency (Novikov & Thorne 
1973). 

We assume that the power available for the production 
of jets can be described as a function of accretion rate and 
spin only: 



Qjct = rj(a)m,c 2 



(2) 



where <5j e t is the jet power and 77(a) is the jet efficiency 
This is a simplification, which ignores the dependence of 
the jet power on the geometry of the accretion flow. For 
more details see MSR11. In these proceedings we only 
show the results using one particular set of jet efficien- 
cies, from the 3D general relativistic magnetohydrody- 
namic simulations of Hawley & Krolik (2006). However, 
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Fig. 1. (Left panel): Best- fit to the local radio LF of the ADAF and QSO components. The black circles represent the 
observed radio LF of LEGs, the empty triangles represent the HEGs. (Right panel): Spin distribution for all SMBHs 
in the loci Universe, obtained by combining the individual spin distributions of ADAFs and QSOs, weighted by their 
space densities (see Equation [5]) . 



in MSR11 we show the results for a set of six different 
efficiencies, and find the results to be robust. 

Finally, to convert from jet power to observed ratio 
luminosity at 151 MHz, we follow the conversion derived 
by Willott et al. (1999): 



Qjct 



= 3 x lcrv 3 



151 



10 28 WHz-isr- 



(3) 



and assume / =20 (see MSR11 for more details, and also 
Cavagnolo et al. 2010). The term / is one of the domi- 
nant sources of uncertainty in our work. We convert from 
151 MHz to 1.4 GHz assuming a power law with index 
-0.75. 

2. Modelling the local radio luminosity function 

Best et al. (in prep.) have classified the radio sources of the 
local radio luminosity function (LF), according to their 
optical spectra. Sources have been classified into high- 
excitation and low-excitation galaxies (HEGs and LEGs, 
respectively), and individual radio LFs have been derived. 

We model these two populations independently. For 
the HEGs, we assume that these represent SMBHs ac- 
creting at a significant fraction of their Eddington limiting 
('QSOs'). The X-ray LF provides us with the space den- 
sity of SMBHs with high accretion rates (Silverman et al. 
2008). The bolometric luminosity, Lboi> can be estimated 
from the X-ray luminosity Lx, via a bolometric correc- 
tion, Cx, so that Lb i—CxLx- Using equations [1] and [51 
the jet power can be estimated using: 



Q 



jot 



rj(a)CxLx 
e(a) ■ 



Hence, given a spin we can infer the jet power from the 
X-ray luminosity. We model the space density of HEGS 
with a given radio luminosity as: 



ISO 



{L v ) = 



dQ 



jet 



dL v 



Pqso(a)(f>(Lx)'^ L da. 



(5) 



The term (/>qso{Lv) is the modelled radio LF of QSOs, 
which we use to explain the HEGs, 4>(Lx) is the X-ray 
LF, and Pqso(o) is the distribution of spins for the QSOs. 
dQjct/dLv is computed from Equation^ For details of the 
derivation of Equation [5] we refer the reader to Section 
4.1 of MSR11. The X-ray LF is integrated above X-ray 
luminosities that can only be reached by SMBHs with 
tti.kJ 10 s M Q . We note that we include a correction factor 
to account for the fraction of QSOs that are optically-thick 
to Compton scattering and hence missed by the hard X- 
ray surveys (e.g. Martmez-Sansigre et al. 2007). 

For the LEGs, we model these as low-Eddington rate 
objects ('ADAFs'). We do not have direct access to an 
ADAF LF. However, the local mass function of SMBHs 
is known (e.g. Graham et al. 2007), so we can use it to 
model the space density of ADAFs. Given a black hole 
mass, to,, the Eddington limiting accretion rate is given 
by TOeC 2 = LeTO./c, where Le is the Eddington limiting 
luminosity, 1.3xl0 31 W M© -1 . Hence, given a black hole 
mass and an Eddington rate, A = to/toe, we can estimate 
the jet power: 



jot = 77(a) Ato e . 



(6) 



Given the space density of SMBHs with a mass to,, 
0(to,), the space density of radio sources powered by 
ADAFs, 0adaf(-M, is given by: 



dQ 



(4) dL, 



jet 



-Padaf(S) 



0ADAF {Lii) = 

P(X)(j)(m,)— — -dTO.do, 
dm' 



(7) 



where Padaf(») is the spin distribution of ADAFs and 
-P(A) is our prior for the distribution of Eddington ratios 
(again see MSR11 for more details). Given our ignorance 
of the distribution of Eddington ratios, we assign a flat 
prior in log space, with uniform probability density in the 
range —8 < log 10 A < —2. 
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Fig. 2. (Left panel): Comparison of the z =1 predicted and observed radio LFs. The solid line represents the observed 
luminosity function (Willott et al. 2001, Smolcic et al. 2009). The dashed line is the predicted QSO contribution, the 
dotted line is the ADAF contribution. (Right panel): The corresponding spin distribution of SMBHs at z —\. 



The only free parameters are the terms describing the 
spin distributions -Pqso(q) and Padaf(q-)- We use the data 
to constrain the best fitting parameters for three different 
spin distributions: a power law, a single gaussian and a 
double gaussian. The bayesian odds ratio is used to choose 
between models. 

Figure [1] (left panel) shows the observed radio LFs 
for HEGs and LEGs, and overlayed are the best-fitting 
4>qso(L v ) and </>adaf(£i/)- The best fitting distribution 
for the QSOs is a single gaussian centred around a=0.00, 
for the ADAFs it is a double gaussian, centred at a=0.06 
and 0.99, and with the high-spin gaussian having an am- 
plitude of 0.79 compared to the low-spin gaussian. 

Hence, we find that the high-accretion rate SMBHs 
(QSOs/HEGs) have typically low spins, while amongst the 
low- accretion rate SMBHs (ADAFs/LEGs) there is a high 
fraction of objects with high spins. 

The typical spin distribution for all SMBHs at z=0 can 
be estimated from the weighted mean: 

-Psmbh(cj) = 

J2 <ftQSo(£i/)-PQSo(a) + J2 </>ADAF(^)-PADAF(q) ^ 

J2 </>qso(-M + S <£adaf(-M 

The right panel of Figure [T] shows -Psmbh(o) in the 
local Universe. 

3. Predicting the z =1 radio luminosity function 

We now test whether our modelling can reproduce the 
radio LF at z =1. From the observed X-ray LF, we can 
model the evolution of the QSOs. 

No such information is available for the ADAFs, how- 
ever. Given that the bulk of the growth of the most 
massive black holes, with m.> 10 8 M Q , ocurred in the 
range 1 JSJ z JSJ 3, make the assumption that there has 
been negligible evolution in the mass function of SMBHs. 
This assumption is supported by the weak evolution ob- 
served amongst radio sources of moderate luminosity (e.g. 



Smolcic et al. 2009). Hence, at z =1 we use the lo- 
cal SMBH mass function, with the same distribution of 
Eddington ratios. 

Figure [2] (left panel) shows the resulting radio LF for 
QSOs (dashed) and ADAFs (dotted), compared to the ob- 
served radio LF (solid line). Without any more free param- 
eters, only the assumption of no evolution for the ADAFs, 
we are able to approximately reproduce the z —1 radio LF 
over six decades of radio luminosity. 

4. The evolution in cosmic spin of SMBHs 

Applying Equation [5] to z =1, we obtain the spin distri- 
bution for all SMBHs at z =1, which is shown in the right 
panel of Figure [2] We see that at high redshift, the fraction 
of SMBHs with high spin is lower than at low redshift. 

It is expected that at z > 1 the mass function of 
SMBHs will begin to decrease significantly in space den- 
sity. This suggests that the space density of ADAFs will 
also decrease, so that we expect the spin distributions at 
higher redshifts to be even more dominated by the a~0 
term. 

The evolution of the spin of SMBHs can therefore be 
understood as a gradual switch between two populations. 
At high redshift, a high-accretion rate but low spin pop- 
ulation dominates (the QSOs). This population decreases 
strongly in space density as z decreases, revealing a second 
population of low-accretion sources, which have a bimodal 
spin distribution (the ADAFs). 

This evolution is best described by the spin distribu- 
tions, as illustrated by the right-hand panels of Figures Q] 
and[H However, for completeness we also show the evolu- 
tion of the mean spin (the expectation value) as a function 
of redshift, in Figure [3] (left panel). This shows a modest 
evolution from (a(z))^0.2 at z ~ 1, to (a(^))~0.35 in the 
local Universe, and when extrapolated to higher redshifts, 
predicts a mean radiative efficiency ~0.065, in excellent 
agreement with observational constraints (e.g. Martmez- 
Sansigre & Taylor 2009). 
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Fig. 3. (Left panel): Evolution of the mean spin inferred from our spin distributions for all SMBHs as a function 
of z. The solid line is the mean, the dashed-dotted line shows the standard deviation (note the -la goes along the 
x axis). (Right panel): Evolution of the mean spin from the cosmological simulation of Fanidakis et al. (2011), for 
m.> 10 8 M . 



5. Discussion 

We have found that the population dominanting the radio 
luminosity function switches from high-accretion rate ob- 
jects with low spins at high redshift, to low-accretion rate 
objects with a bimodal distribution of spins, with approx- 
imately half the population having maximal spin. 

This can be explained by the effects on the spin of 
SMBHs by the two mechanisms by which these grow: ac- 
cretion and mergers. Continous accretion along one plane 
would lead to SMBHs having essentially maximal spin. 
However, it has been recently suggested that SMBHs do 
not accrete in this way, but are rather subject to 'chaotic 
accretion', where matter falls in from different directions. 
Approximately half of the matter has angular momentum 
in the opposite direction to the SMBH, and the final spin 
is therefore, on average, close to 0. Hence accretion pro- 
vides a mechanism for spinning SMBHs down (King et al. 
2008). 

Major mergers of SMBHs, on the other hand, are likely 
to increase the final spin (e.g. Rezzolla et al. 2008). When 
two black holes of similar mass merge, the angular mo- 
mentum of the final orbit will be significant compared to 
the final mass, so that the spin of the coalesced black hole 
will be large. 

The evolution of the spin distributions inferred in our 
work is in good agreement with a picture where at high 
redshift the spins are typically low due to accretion. Major 
mergers will occur both at high and low redshift, but at 
high redshift there is a larger supply of cold gas available 
(Obreschkow & Rawlings 2009), so that accretion will spin 
the SMBHs back down. 

At low redshift, the cold gas is running out, so after 
the mergers the breaking mechanism is no longer present. 
If a large fraction of the SMBHs have undergone a recent 
major merger, then a large fraction of the SMBHs will 
retain a high spin. 

A detailed cosmological simulation of the growth of 
SMBHs was performed by Fanidakis et al. (2011), which 



included the effects of black hole mergers as well as chaotic 
accretion. The left panel of Figure [3] shows the mean spin 
from their simulation, as a function of redshift. The simu- 
lation shows an almost identical evolution to that inferred 
from our work, and a very similar, large variance. 
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